Loss of coherence and dressing in QED 
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The dynamics of a free charged particle, initially described by a coherent wave packet, interacting 
with an environment, i.e. the electromagnetic field characterized by a temperature T, is studied. 
Using the dipole approximation the exact expressions for the evolution of the reduced density ma- 
trix both in momentum and configuration space and the vacuum and the thermal contribution to 
decoherence, are obtained. The time behaviour of the coherence lengths in the two representations 
are given. Through the analysis of the dynamic of the field structure associated to the particle the 
vacuum contribution is shown to be linked to the birth of correlations between the single momentum 
components of the particle wave packet and the virtual photons of the dressing cloud. 

PACS numbers: 03.65.Yz, 03.70.+k, 12.20.Ds 



I. INTRODUCTION 

Decoherence consists in the destruction of coherences 
present in the initial state of a quantum system due to 
the interaction with external degrees of freedom [l| . De- 
coherence is associated to the increase of entropy and the 
loss of purity of the initial state of the system [2] • 

The environment may be regarded as monitoring cer- 
tain properties of the quantum system through the in- 
teraction with the system itself Q. Not all initial quan- 
tum states are equally fragile to this interaction: often 
there are relatively robust states with respect to it, called 
"pointer states" |4|. Experimental evidence of this envi- 
ronment induced decoherence has also been recently re- 
ported [l! [Ell- 
in the case of a particle, either free or in a potential, 
linearly coupled to the environment modelled as a bath 
of harmonic oscillators at temperature T, several stud- 
ies of decoherence processes have already been reported 
US [HE! EE In thcsc studies both the Hamilto- 

nian approach and functional techniques have been used. 
It has been shown that, starting with the particle and the 
bath described by a factorized density matrix, it is possi- 
ble to distinguish two characteristic contributions to the 
decoherence: the first related to the thermal properties of 
the bath and the second, independent of temperature, to 
the zero point fluctuations of the oscillators of the bath 
EB|- Decoherence has been shown for charged particles 
initially described by a wave function made of a coherent 
superposition of two moving wave packets to be linked to 
the emission of Bremsstrahlung E3 • 

Here we want to investigate the role played by radi- 
ation emission and entanglement with field degrees of 
freedom, on the decoherence induced on a free charged 
particle by its interaction with the electromagnetic field 
at temperature T that plays the role of environment. 
In particular we shall study the decoherence among the 
components of an initially Gaussian free wave packet rep- 
resenting the particle by analyzing the evolution of the off 
diagonal elements of the particle reduced density matrix. 



We shall start, as is typically done [18|, [l9|, [2fJ, [2l| , from 
decoupled initial conditions which correspond to the ab- 
sence of initial correlations between the system and the 
environment. Using physical approximations, we'll re- 
duce the particle-field interaction to a simple analytically 
solvable model. 

We will focus mainly on two aspects. The first one 
is the analysis of the dynamics in two different basis. 
The aim is of evidencing clearly how the change of rep- 
resentation gives place to different relative importance 
of various effects induced by the coupling with the en- 
vironment. The second one is the study of build up of 
quantum correlations between the system and the en- 
vironment. Using the fact that our model Hamiltonian 
allows exact treatment it is easy to show in detail the 
mechanism linked to the part of decoherence indepen- 
dent of temperature. This will be done by investigating 
the time dependence of the effects on the particle due to 
vacuum fluctuations, such as the dressing, and through 
the analysis of bath dynamics without the use of approx- 
imations as the Markovian one. 

The paper is organized as follows. In Sec. [Til we de- 
scribe the approximations adopted that transform the 
Hamiltonian into a linear form amenable to exact treat- 
ment. In Sec. IIIII the particle density matrix is obtained 
both in momentum and real space. In Sec. IIVI we ana- 
lyze the dynamics of the field structure, evidencing the 
relationship between the part of decoherence induced by 
vacuum and dressing process. In Sec. [V] we summarize 
and discuss our results. In Appendixes 1X1 [Bl [Cl we have 
collected most of the calculations to make more readable 
the main body of the text. 



II. MODEL 

The system under investigation is a free spinless par- 
ticle of mass mo and charge e moving at initial velocity 
Vq, interacting with the electromagnetic field in thermal 
equilibrium at temperature T. The particle is initially 
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described by a coherent wave packet, whose initial width 
is assumed to be small with respect to the relevant wave 
lengths of the electromagnetic field. The interaction be- 
tween the system and its environment is described by 
the non relativistic minimal coupling Hamiltonian with 
an upper cut off frequency fl corresponding to a wave 
length such that the dipole approximation may be ap- 
plied dl. 

The adoption of dipole approximation is standard in 
the treatment of free particle decoherence fill . fl2l . [l3| 
but it limits the validity range to times of the order of 
To = c/voil, where c is the light speed and vq is the ini- 
tial velocity of the particle. This limitation can be made 
less strong by using a "moving dipole" approximation 
which consists in substituting the particle position op- 
erator r by a parameter r t indicating the average wave 
packet position at time t. In absence of interaction this 
is given by rt = ra + vot, ro being the initial position of 
the particle. It is possible to check the consistency of our 
choice by comparing r t — r with the particle average dis- 
placement in presence of the interaction, (q)t, given by 
Eq. (|41|) . In fact their difference is smaller than the wave 
packet width for times less than the ones where moving 
dipole approximation can be applied (sec Sec. [Vj) . Our 
results are valid until a time such that because of the 
spreading the wave packet width becomes of the order of 
the minimal wave length involved in the treatment [l2l |. 
The contribution to the spreading of the wave packet due 
to the interaction can be shown to be for small value of a 
(see Eq. (Unj)) negligible with respect to the free evolution 
for small times. Taking an initial wave packet of mini- 
mum indetcrmination, using Eq. (|46[) for the free spread- 
ing and the dipole approximation condition Ar <C c/Cl, 
we get Td « £7 _1 (moc)/Ap with t c i ^> tq. 

The potential vector in the Coulomb gauge is given by 

M 



A(r)=5> fe 

k,3 



I 2ttHc 2 / f 



(1) 



where e k j are the polarization vectors (j = 1,2) of the 
mode k of frequency u>k, periodic boundary conditions 
are taken on a volume V, h is the reduced Planck con- 
stant and and & k j are the annihilation and creation 
operators of the field modes satisfying the commutation 
rules [&k.j, ji] = fikk'Sjj'- The non relativistic minimal 
coupling Hamiltonian in the " moving" dipole approxima- 
tion is 



where we have used p = ^pP&p, <J V = \p){p\ is the pro- 
jection operator on the momentum p and the potential 
vector is calculated in r t . 

In Eq. ((5]) the term quadratic in A, which is physically 
linked to the average vibrational kinetic energy due to 
vacuum fluctuations p3j . can be exactly eliminated by 
a canonical transformation of the Bogoliubov Tiablikov 
form [JJj]. Here it will be simply neglected because it 
can be shown as usual to be very small compared to the 
linear term. 

Thus, using Eqs. {TJ) and the Hamiltonian reduces 
to the form 



# = E ^ b v + E fe"* a L a w 



afcje 



kr t 



p,k,j 



) , (3) 



with the coupling coefficients given by 



„ e / 2ttH 



(4) 



Here, in contrast to other phenomenological models [24(, 
the coupling coefficients and the spectral field proper- 
ties are assigned, which allows to analyze the dependence 
of the decoherence development on physical parameters 
such as the mass and the charge of the particle. 

The Hamiltonian of Eq. describing the interaction 
between the system (particle) and environment (electro- 
magnetic field) is now treated exactly. 



A. System evolution 

In the interaction picture, introducing the time order- 
ing operator the unitary time evolution operator is 



U(t) = T^ exp 



fiJo 



dsHi(s) 



(5) 



where, from Eq. ((3J, the interaction Hamiltonian at time 
t is given by 

+a fcj e-' i( ^- fe -^V fc - r °) . (6) 



H 



1 



eA(r t ) 



p 



e 



2mo 

E£-^ - • E^ fea ^ a ^ 



E^4,a,j 

o 

A (r t ) 



^a p p ■ A{r t ) . 



2moc 



The commutator of the interaction Hamiltonian at two 
different times is equal to 

[£}(«),#/(«')] = -2i£ & P sm[(uj k - k ■ v )(s - s')] 

P,k,j 

, 9 . (7) 
W where we have used a p & p i = <r p 5 PP '. Because the com- 
mutator ([7]) commutes with the interaction Hamiltonian, 
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it is possible to give an exact expression for the evolu- 
tion operator [3, Il7l| using the Cambell-Baker-Hausdorf 
formula: 



U(t) = cxp 



x exp 



Jo d ' Jo i "' [H ' {s) ' "'Ml* -""> 



■ exp 



Yl £(p> f ) a p 



1 

L P 



x exp 



at ,a fc - a fe: 



(8) 



where 



1 _ p i{u k -k-v )t 

a k = - e -ife-T-o 

h(oj k -k-v ) 



(9) 



The term £(p,t), present in the above phase factor, is a 
number depending on the momentum p and on the time 
t, as it is shown in Appendix |A"1 



III. REDUCED DENSITY MATRIX ANALYSIS 

The analysis of the decoherence of an initial coher- 
ent wave packet will be conducted by examining the be- 
haviour of the reduced density matrix elements. 

As initial condition we take a state with no correla- 
tion between the particle and the electromagnetic field. 
To this condition corresponds a decoupled initial density 
matrix of the form 



p(0) = p S (0) ® PF ■ 



(10) 



where ps(0) represents the initially coherent wave packet, 
while the field is taken in a thermal state at temperature 
T described by pp = cxp(—f3Hp)/Zp, with (3 = 1/ksT, 
ks the Boltzmann constant, Hp is the Hamiltonian of 
the field and Zp the field partition function. 

In Eq. ^ the projection operator a p commutes with 
.ff, thus the particle's momentum is a constant of motion. 
This implies that momentum space provides a robust ba- 
sis that allows to investigate easily the decoherence de- 
velopment. Successively we shall consider the coordinate 
space to see how the loss of coherence shows up in real 
space. 



A. Momentum space 

In the momentum representation the ini- 
tial particle density matrix becomes Ag(0) = 



E^/f (0)b(0))(p'(0)|. 
arc given by 



Its elements at time t 



P P s P (t) = (p(t)\ps(W(t)) (11) 
= ( P (t)\tr F {U(t)ps(0) ® ppU-^tWit)) , 

where \p(t)} is an eigenstate of the momentum operator 
at time t. 

Indicating with \4>(t)} an arbitrary field state we obtain 
U(t)(\ P (t)) ® |0(t))) = \p(t))exp[it(p,t)] 

® ex PXXj (4j°* " ^ a k) w*)> . ( 12 ) 

where use has been made of the fact that the application 

of the operator exp i^2 p £,(p,t)a p of U(t) on the state 

\p(t)) <& |0(O) leads to the factor exp [i£(j>, t)]. This factor 
doesn't depend on the environment state but only on the 
associated momentum. 

We have already seen that with the Hamiltonian Q 
the particle momentum is a constant of motion. The 
states \p(t)) are stationary with respect to the interac- 
tion and different momenta can't be connected by the 
time evolution operator. Then, in Eq. (|11[) . in the mo- 
mentum representation form of ps (0), only the term 

f&iPMowm = P p /m P (t)){p'(t)\ expHv - 

p' 2 )/2mofi] contributes to the reduced density matrix 
evolution. Thus, Eq. (jlip can be written as 



P P s P (t) 
tip 



exp < 1 



t;(p,t)-{;(p',t) 



t{p 2 -p' 2 ) 



2mnh 



[ k,j 



IjOLk ~ &k,jU k pp 



(13) 



where we have used the property of ciclity of the trace. 
We can rewrite this last expression as 



P P s P (<) = P P s P (0) exp -T** (t) + (t) 



(14) 



where we have introduced the decoherence function, typ- 
ically used in literature [H, [13] , as 

r**'(t) = -Y, lnt1 ^ { cx p ( a L<f - h,^'*) pf} 



(15) 



with = (gl tj - glj)a k , and the function (t) 



$ p ' p (t)=Z( P ,t)-Z(p',t) 



t(p 2 -p' 2 ) 
2moH 



(16) 



that includes the phase term t) — £,(p', t) and the free 
evolution term. 

The decoherence function describes in a direct way 
the appearance of decoherence. In fact, the increase of 
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pp,p (£) for p =/= p' gives rise to a decrease of the off di- 
agonal elements of the reduced density matrix, that is it 
leads to the destruction of coherences among the differ- 
ent momenta in the initial wave packet. Moreover, the 
expression of 7^'? shows that at p = p ' the decoherence 
function is zero and then that the populations are con- 
stant in time. This may be expected because, as shown, 
the dipolc approximation leads to momentum conserva- 
tion. 

For our model the calculation of the explicit form of 
the decoherence function T p,p (t) and the phase factor 
t) — t) is reported in Appendix lAl 

Eq. (|A20[) shows that the decoherence function T pp (t) 
increases quadratically with the vector difference of the 
momenta p — p' . Therefore there is decoherence in the 
off diagonal elements also within the same energy shell. 
Introducing the spectral density, 



J(u) = 



2a (p-p') 



3- 



-,2 n 2 



■ lo exp 



(-5) 



(17) 



containing the frequency dependent part of Eq. (|A20[) de- 
riving from the coupling coefficients and the density of 
the modes at frequency oj, T p ' p (t) can be rewritten as 



dw J(u>) 



(1 — cos Ult) 



coth 



2k B T 



(18) 

Below the cut off frequency ft, J(to) depends linearly on 
lu, this is typical of an Ohmic spectral density which gives 
rise to frequency-independent damping [l7j | . This damp- 
ing gives rise to a loss of coherence between different mo- 
mentum cigenstatcs but not to dissipation, which is ab- 
sent because the interaction Hamiltonian commutes with 
the momentum operator. 

In Appendix [3] it is shown that it is possible to sep- 
arate in the decoherence function T p ' p (t) the effects of 
vacuum fluctuations, T p ' p c (t), and of thermal contribu- 
tion, T% p '(t), as T p ' p '(t) = T p £ (t) + T p tk p ' (t) . Extracting 
the dependence on the momenta we rewrite the decoher- 
ence function as 



T p - P '(t) = T(t)(p-p'Y 



(19) 



2a 



In 



Vi + ft 2 * 2 



sinh(i / rp ) 
t7rF 



3tt 



2 9 



(p-p'Y 



with T(t) the decoherence factor and a = e 2 /he a dimcn- 
sionless coupling constant. For the two contributions we 
write 



T p v *(t)=T vac (t)(p-p' 



2a ln\/l + ft 2 i 2 
37r unkc? 



ip-p'f 



(20) 



with T vac (t) the vacuum decoherence factor and 



T™{t) = T th {t){p-p'f = 



2a 



In 



ih(t/r F ) 



t/T F 



3tt 



,2„2 



(P-P 



'\2 



(21) 



with T t h(t) the thermal decoherence factor and rp = 
h/nksT a characteristic thermal time. The expression 

for Y P £ (t) is obtained under the condition fc#T <IC Ml. 
If Ml m 10 _2 m c c 2 , m c indicating the mass of an electron, 
the above condition is well verified at ordinary conditions 
(T < 10 7 ^). 

Eq. (UHJ) shows that T p p (t) increases faster with time, 
the difference p — p' and the coupling constant a. 
From Eq. (|A12|) we obtain: 

Ofy f, 2 — f,' 2 

£(p, t) - Z(p', t) = — ^—s^-(SU - arctanftt) , (22) 

which depends only on the energy difference between the 
components of momentum p rather than on their vector 
difference. Separating the dependence from momenta we 
introduce from Eqs. (fTTJ)) and (12"2")) the global phase factor 
$(t) as 



<I> P ' P (t) 



„'2\ 



2a (fit — arctan fit) 



t 



37TTOqC 2 



2mofi 



(23) 
(p 2 -p' 2 ) ■ 



We observe that $(i) doesn't depend on the initial state 
of the field and that in absence of interaction it represents 
the phase free evolution, <f>(i) = —t/2mH. 

Using Eq.(P|) for T p > p ' (t) and Eq. (J23) for $ p > p '(t), 
we can rewrite the particle density matrix elements of 
Eq. (HP as 

p p s p \t) = p^'(0) exp[-r(i) (p - p') 2 +l $(t) (p 2 p' 2 ) . 

(24) 

To discuss the time evolution of the reduced momen- 
tum density matrix elements it is useful to use simplified 
expression for T(t) and <f>(t) for different times easily ob- 
tainable from Eqs. ([19]) and l[23l: 



r(t) « ^ 



T th (t) 



2a fl 2 t 2 
272 — o~ 



2a 



in ftt, 



2a t 



ft" 1 < i < rp , 

i > rp , 
(25) 



and 



S(f) 



3 + 3 



2 a flH 
Zitm^c 2 3 



2nif)h 



2a n t 



37TTOqC 2 



i< ft" 1 



t > ft" 



(26) 



We observe that the form of the decoherence factor T(t) 
leads to a time behaviour for the reduced density matrix 
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elements analogous to the one obtained for an ensemble 
of two level systems linearly interacting with a bath of 
harmonic oscillators Q. In our case there is an explicit 
expression of the coefficients in terms of the parameters 
of our system. 

Using in Eq. (f24|) the approximated expressions of T(t), 
in the three time zones of Eq. (|25[) . and the expansion 
c~ x ?s 1 — x for x <C 1 we obtain 



// (*) 



1 - 



2a (p-p'f n 2 t 2 



3tt 



,2„2 



t < O" 



cxp 



(P-P') 2 



2a (p-p') 2 * 

37T ?71qC 2 Tp 



t > T F . 

(27) 



Eq. (|27[) shows that the off diagonal elements of p P g P (t) 
evolve from the initial value for small times with a 
quadratic trend, for intermediate time with an hyper- 
bolic and at large times with an exponential one with 
the rate 2a(p — p') 2 /SnrriQC 2 . 



1. Vacuum and thermal contribution: decoherence times 

It is possible to use the approximated expression of 
Eq. (|25p for T(t) to evidence the time regions in which 
vacuum and thermal contribution dominate. It comes 
out that the vacuum contribution prevails for t <C Tp 
while the thermal contribution dominates for t > Tj?. 
The transition time, r p , at which the two contributions 
are equal can be found imposing In Qt p = T p /rp. This 
time doesn't depend on ( = 2a(p — p') 2 /SirrriQC 2 . For 
example for Q ps lO 1 ^ 1 (Ml « m o c 2 /100) and T = IK 



we have Tp ps 2.34 • 10 



from which we find r„ 



In Fig. 1 the behaviour in time of T%£(t) and T p th p ' (t) 
is shown as a function of physical parameters present in £. 
It shows that if T p ' p c (t p ) > 1 then vacuum contributes ef- 
fectively to decoherence, otherwise only the thermal con- 
tribution will be effective. 

In the range where the vacuum contribution dominates 
(t Tp) there are two different typical time dependen- 
cies. In the first one (t <C f2 -1 ) the increase of decoher- 
ence is fast while in the second one (t fiP 1 ) it slows 
into a logarithmic dependence. Fig. 2 represents the 
time development of exp[— TPf c (t)] as a function of the 
coupling constant a, showing that by increasing a and 
fixed p — p', we observe a decay of matrix elements due 
to the vacuum contribution faster in time. 

We distinguish two different characteristic times of the 
decoherence process relative to the vacuum 



1 

7~vac — GXp 



2a (p — p ' 



(28) 



C=0.25 



f — i 



P,r- 
F — i 




C=o.oi 



50rj 



FIG. 1: Figure shows T v v 'g c (t) and (t) as a function of 
C = 2a(p - p ') 2 /3-Kmlc 2 , in the case T = 3WK. 




FIG. 2: In figure it is plotted the time development of 
exp[— Ty'ac (t)] (|20[) as a function of the coupling constant 
a = e 2 /he in the case \p — p'\/moc = 0.1. The time is taken 
in unit of tQ and the range of a is chosen to visualize the 
vacuum effect on the decoherence. 



and to the thermal contribution 



Tth = Tf 



3tt m 2 c 2 
2a (p - p'f 



(29) 



These characteristic times have the same form of those 
obtained for the decoherence of the interference pattern 
in [13. 

The mass and charge parameters mo and e, appearing 
in T vac and Tth, are arbitrary. The only restriction is that 
they refer to a body that can be treated as a point like 
particle within the dipole approximation. For example, 
these parameters could represent the mass and the charge 
of a highly charged nucleus or even of a macroscopic body 
of linear dimensions small enough, and therefore a is a 
free parameter. 

Let's observe that the time at which vacuum and ther- 
mal decoherence are effective, depending of the value of 
the coupling constant a, fall inside the time of validity 
of our model. 
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2. Analysis of Ap(t) and l p (t) 

The above results are independent from the structure 
of the initial reduced density matrix elements p P g P (0). 
Now we specialize these results to the case of an initial 
Gaussian wave packet of spatial width Ar 



fff (0) = N exp ■ 



3[(p-Pq) 2 + (p'-Pq) 2 ] 
AAp 2 

ro ■ (P-P r ) 
~ l ft 



(30) 



with Ap the width in the momentum space. p the initial 
average momentum of the particle, N = (v / 3/v / 27rAp) 3 
the normalization factor and ArAp = 3ft/2. 
Substituting the gaussian wave packet of Eq. (|30|) in the 
reduced density matrix at time t of Eq. (pM|) , this can be 
put under the form 



fff (t) =Nexp 



x exp 



- r(*) + 



2Ap 2 
3 

8Ap 



(31) 



{P-P 



I \2 



3(p + p Q 
8Ap 2 



>Y2 



+i$(t)(p 2 - p' 2 ) + 



3p '(P + P') .^o-(p-p') 



2A P 2 



A way to quantify the degree of loss of coherence of the 
wave packet is through the coherence length l p (t) [l|, de- 
fined as the width of ps(t) along the main skew diagonal, 
meaning the region inside which the coherence between 
momenta has not been yet destructed at time t. l p (t) 
may be compared with the width of ps(t) along the di- 
agonal that measures the wave packet width at a time 
Ap{t), given by 



Ap(t) 



< p 2 > t - < p >\ = Ap, 



(32) 



where we have used < p 2 >= tr(ps{t)p 2 ) = Ap 2 +p^ and 
< p > t = tr{ps{t)p} = p Q . Because Ap(t) is constant 
the wave packet doesn't spread with time in momentum 
space. 

The coherence length l p (t), proportional to the inverse 
of the square root of the coefficient of {p—p ') 2 in Eq. (|3"Tj) 
[1, is: 



Ap(t) 



y/1 + 8A P 2 T(t)/3 



(33) 



To quantify the effective loss of coherence in the wave 
packet we study the ratio l p (t)/Ap(t). This quantity 
gives a measure of the relative width of the reduced den- 
sity matrix off the diagonal compared with the width 
along the diagonal. Using Eq. (|3"3"|) and the explicit form 
of r(t) for t > Tp, given by Eq. (|2l)|) . we obtain 



l p (t) 3moc 



Ap(t) AAp 



TTTp 

at 



(34) 



where t* = (3moc) 2 TTTF / (4:Ap) 2 a. Being Ap(t) constant 
in time, Eq. (|3~4"]) shows that the coherence length for large 
times decreases going to as l/y/t for t — ► oo. The deco- 
herence process in momentum space is thus characterized 
by a complete decay of the off diagonal elements of the 
particle density matrix for large times while the popula- 
tions remain constant. 

This kind of behaviour of the reduced density matrix 
is shown in Fig. 3 obtained from Eq. (|3 1|) . 




0. 2 



FIG. 3: In figure it is represented the absolute value of the 
normalized density matrix Z = \pg" a (t)\/N x in one dimen- 
sion, N x = l/y/2TrAp x , with po = and Ap x /moc ~ 0.1. On 
the top it is t = while on the bottom it is t = 3r va c, where 
T vac is taken for p x — p' x = Ap x . 



B. Coordinate space 

Our analysis is now extended to real space in order to 
describe spatial decohcrencc in more complex situations 
such as Young interference or Schrodinger cat states se- 
tups. We expect that changing representation the dy- 
namics induced by the interaction with the electromag- 
netic field will appear more complex than in momentum 
space. As shown, in fact, it provides a basis of pointer 
states which allows a simple analysis of the process. To 
investigate the effects in the real space we need the re- 
duced density matrix in the configuration space. It can 
be obtained from the corresponding momentum space 
reduced density matrix by performing a double Fourier 
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transform: 



d'pdypY (i)exp 



(p ■ r — p' ■ r') 



(35) 



Taking the Gaussian wave packet described by p P g P (0) of 
Eq. (|30| . the transform can be explicitly performed and is 
given in Appendix[Bl The spatial reduced density matrix 
p T g T (t), given by Eq. (|B8j) . can be rewritten as 



P7 (t) 



x cxp 



x cxp 



NAp 3 



cxp 



.Ar 2 + 6r(t)fi 2 Po-(g-Q') 



Ar(i) 2 



3[(g-(g) t ) 2 + ( 9 '-(g) t ) 5 



Ap' 



4Ar(t) 2 

2 [-r(t)(g-q') 2 - ^)(? 2 -9' 2 )]} 



(36) 



Ar(t) 2 

where: q = r — To is the displacement from the initial 
position and its average at the time (q)t is given by 

(q)t = Hp s {t)q} = -2p ^(t)H , (37) 

$(t) and r(t) are defined by Eqs. 03]) and {23j) while 
Ar(t) is the spatial width of the wave packet at time t. 
From Eq. |[B8|) we get 

(g 2 ) t = tr{p s (t)g 2 } (38) 

9<J>(t) 2 fi 4 



Ar' 



6T{t)H 2 



2fc2 



Ar 2 



and thus, using also Eq. ([57]) . Ar(i) is given by: 
Ar(i) = Aq(t) -- 



An /l 



6r(t)fi 2 9$(t) 2 fi 4 



Ar 



Ar 



(39) 



Ar(i) is at t = 0, Ar(0) = Ar, and increases with time. 
p r ^ r (t) can be obtained by its form at t = 

.p a -(q-q') 3(q 2 + q /2 Y 



r y NAp 3 
Ps (°) = ex P 



4Ar : 



(40) 

replacing the initial width of the wave packet Ar with 
its value at time t, Ar(t), multiplying by (Ar 2 + 
6T(t)h 2 ) I Ar(t) 2 the phase factor in the first exponent, 
centering the wave packet in the average displacement 
(q}t in the second exponent and multiplying by an expo- 
nential factor which gives an increase of decohcrence and 
a phase variation analogous to the factor appearing in the 
reduced momentum density matrix elements of Eq. (|24|) . 



1. Time dependent dressing 

The average of the operator q at time t, given by 
Eq. ([37]) and using the explicit form of $(i) ([23]), is 



Pot 



4ah _ arctannn 

37rmoC 2 



(41) 



From this equation the average velocity of the wave 
packet is 



df 



[q)t 



m 



1 



Snmoc 2 1 + n 2 t 2 



(42) 



As observed before, p is a constant of motion, instead the 
velocity q = [p — eA(ro)/c]/mo it is not because it does 
not commute with the Hamiltonian ([3]). This may be 
related to the fact that, starting with uncoupled initial 
conditions, the charged particle is subject to time depen- 
dent dressing by the transverse photons. This increases 
its mass while p remains constant. The mass variation 
can be obtained casting Eq. (|42[) in the form 



(h)t 



m(t) 



Po_ 

too 



1 



Sm(t) 



to 



(43) 



where m(t) = Too + 5m(t) is the mass at time t being the 
mass increase 5m(t) given by 



6m(t) = 



2+2 



AahVL fl 2 t 

3ttc 2 i + n 2 t 2 



AaHn 
3irc 2 



2+2 



AahVL 
3tt c 2 



t > n~ 



,_. (44) 

For t <ti £1 1 5m(t) increases quadratically [26j while for 
t ^> f2 _1 coincides with the usual total mass variation 
due to the interaction with the electromagnetic field [13] • 
We observe that the equation of motion ([41"]) . from 
which wc derived the expression for the mass increase, 
is related only to the total phase factor <!>(£) and is then 
temperature independent at first order in Vq / c. 



2. Analysis of Ar(t) and l r (t) 

The mass variation due to dressing is relevant if one 
wishes to compare the evolution of the wave packet width 
in the absence of interaction, Ar(t) , with its expression, 
Ar(i), in the presence of interaction. In the last case we 
have from Eq. 



ApH 2 4fi 2 $(i) 2 6T(t)H 2 
Ar(t) = An, 1 + + (45) 



A; 



with T(t) and defined by Eqs. ([19]) and ([23]). Putting 
e = we obtain the well known expression for the free 
spread [25[ 



Ar(t)° = ArJl 



Ap 2 i 2 1 



Ar 2 TOq 



(46) 



(q)t, given by Eq. ([37]) . can also be obtained by integrat- 
ing Eq. (22) 



f* 1 

(q)t = -2p fim) = p / 77a ' 
Jo m \ l ) 



(47) 
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Thus we can identify 



-2m(t) 



m(t') 



dt' 



m(t') 



(48) 



rn(t') 



being the time average of l/m(t') over the time 

t. The width of the wave packet at time t (|45[) can be 
thus rewritten as 



Ar(t) 




1 



m(t') 



6T(t)H 2 
Ar 2 



(49) 



Eq. (|49p shows that, starting from uncoupled condition, 
the interaction with the electromagnetic field induces dif- 
ferences with respect to the free evolution Ar(t)°. The 
first one consists in the replacement of the inverse of the 
initial mass 1/mo by (l/m(t')) t and may be attributed 
to the ^-dependent dressing. This effect is due to the vac- 
uum fluctuations and is related to the total phase factor 
<!>(£), the mass increase leading to a rate decrease of the 
width with respect to the free case. The second effect is 
given by the term within the square root 



6T(t)h 2 
Ar 2 



1 Aah 2 



In 



y/i + n 2 t 2 



sinh(i/TF) 



t/T F 



(50) 

It always leads to an additional increase of the width of 
the wave packet. It contains both the effect of vacuum, 
represented by the term \/l + £l 2 t 2 , and of the thermal 
field represented by the term sinh(£ / Vp) / (t / Vp) > being 
this last term for T = equal to 1 (tf = fi/7r/cgT). 

The comparison of the amplitudes of the vacuum and 
thermal terms in time may be obtained using the forms 
of the coefficients T(t) and <f>(t) given by Eqs. (f25|) and 
([26]) for small (t <C fi _1 ) and large (t 3> tf) times . 
For small times the total effect is that the width of the 
wave packet results larger than in the free case. For large 
times, instead, the additional term becomes negligible 
and the spreading is slower than in the free case because 
the increasing of mass. 

The space coherence length l r (t) represents the typical 
distance for which it is possible to have constructive inter- 
ference among different parts within the wave packet. It 
can be read directly from the coefficient of {q — q 1 ) 2 term 
of the reduced density matrix written under the form of 
Eq. (|B9[) , being in fact proportional to the inverse of this 
coefficient p|: 



l r (t) 



Ar(t) 



v /l + 6fi 2 r(i)/Ar 2 



(51) 



Using Eq. ([39]) for Ar(t) it results that l r (t) increases 
with time, while, analogously to what happens in mo- 
mentum space ([3"3"j) . l r (t)/Ar(t) decreases with time be- 
cause T(t) increases with time ([T9|). In absence of inter- 
action T(t) is equal to zero and therefore the free space 
coherence length, l r (t)°, is always equal to the width of 
the wave packet which increases coherently in time due to 



the well known free spread ([46)) . The coupling with the 
field induces an evolution of l r (t) different from Ar(t). 
Using Eqs. (J33J) and ([51]) and ArAp = 3H/2, it follows 
that l r {t)/Ar(t) = l p (t)/Ap(t) and therefore Eq. ^ de- 
scribes also in the coordinate space the behaviour of the 
coherence length with respect to the width of the wave 
packet for large times. This equation shows that the ra- 
tio decreases to zero as 1/Vt for t — > oo describing a loss 
of coherence also in the configuration space. 

Another interesting aspect to investigate is the be- 
haviour of l r (t) with respect to its evolution in the free 
case l r (t)°. Using Eq. gBJ) for l r (t)° and Eqs. (J3TJ> and 
(|49|) we can put the coherence length in the form: 



l 2 (t) = Ar 2 + Ap 2 t 2 



1 



m(t')/ t Ar 2 + 6T(t)ff 



<l 2 (t)°. 



(52) 



Eq. (l52|) shows that dressing induces a slower increase of 
coherence length due to the mass increase, but always 
maintaining the coherence, while vacuum and thermal 
field induce a destruction of coherence in space such that 
the coherence length is lower than in the free evolution 
case. 

In the momentum space we obtained a simple dynam- 
ics: the width of the wave packet remains constant while 
the coherence length decreases with respect to its initial 
value going to zero. In coordinate space, instead, dif- 
ferent factors contribute to the dynamics: free evolution 
contributes to the coherent increase of the width of the 
wave packet coherently and therefore of the coherence 
length; the particle time dependent dressing of the parti- 
cle slows this increase; finally vacuum and thermal field 
induce a loss of space coherence such that the value of the 
space coherence length in presence of the interaction is 
always lower than its value in absence of the interaction. 



3. Linear entropy 

The dynamics of our system is described by the re- 
duced density matrix time evolution as a transformation 
from the pure initial state (jTUJ) into a statistical mixture 
([2"4"| . The time dependence of this process, that implies 
a loss of information on the system, may be described 
by the so-called linear entropy, Sn n [131 • ^ has been an- 
alyzed in the case of localization by scattering, to mea- 
sure how strongly the environment destroys coherence 
between positions by dclocalizing phases, finding a linear 
departure in time from the initial value describing a 
pure state [l[. Using its definition we obtain here 



Sim = tr(p s - p%) = 1 



l + 6T(t)h 2 /Ar 2 



(53) 



which describes the loss of purity of the initial state. It 
is interesting to note that in the case of initial Gaussian 
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wave packet [271 ] . Sim is directly connected to a dimcn- 
sionless measurement of the decoherence given by the ra- 
tio between the decoherence length and the wave packet 
width. This ratio coincides both in the p and r represen- 
tations (j33|) and (|5l"|) and using Eq. ([53]) may be expressed 
as 



Si 



lin 



lp(t) _ l r (t) 



Ap(t) 



Ar(t) 



(54) 



Using Eq. (|33[) and the approximated form of T(t) for 
small times given by Eq. (|25p , we find that Su„ at the 
beginning evolves quadratically from the initial value 
corresponding to a pure state, then slows and finally (|3~4"|) 
goes to 1 for t — > oo as 1 — 1/ 

In Fig. 4 the time development of the linear entropy 
is plotted as a function of the coupling constat a. The 
figure shows clearly that the increase of this quantity 
towards 1 depends strongly on a, that is on the charge 
of the particle considered. 




1/137 



FIG. 4: In figure it is shown the behaviour of Su n in time as 
a function of a. 



IV. INTERPRETATION OF VACUUM 
INDUCED DECOHERENCE 

The temperature independent part of decoherence is 
represented by T^f c (t) of Eq. ([2U)l . In the following we 
shall analyze the processes that contribute to (t) . 

In the case of a charged particle, initially described 
by a wave function made of a coherent superposition of 
two moving wave packets, it has been previously shown 
fl6| that Bremsstrahlung radiation induces decoherence 
decreasing the visibility of the interference pattern that 
results from their overlapping. The reason is that in its 
trajectory the particle is subject to a sudden change of 
the 4-momentum and in this process it is radiated as 
Bremsstrahlung photons the energy [28j ] : 



where I(v, v') is a function of the initial and final ve- 
locity and k max is the wave vector corresponding to the 
frequency equals to the reciprocal of the time scattering 
during which the 4-momentum changes. This energy re- 
sults, as also the decoherence function, proportional to 
a. Thus Bremsstrahlung may be hold responsible of de- 
coherence. 

In our system the particle is also subject to a change 
of velocity during the dressing process with the emission 
of Bremsstrahlung photons. These could be held respon- 
sible of the temperature independent loss of coherence 
between the momentum components of the wave packet. 
However the radiation energy emitted in the unity of time 
from the accelerated charged particle during the dressing 
process can be estimated as [29j : 



E B oc 



e 2 (q) 2 



(56) 



with (q) being the average acceleration of the particle. To 
obtain (q) during the dressing we take the time derivative 
of Eq. (01 : 



p AahVL 2n 2 t 
mf 3ttc 2 (T+nap 



(57) 



(55) 



Substituting this last equation in Eq. (|56[) the estimated 
emitted energy per unity of time results proportional to 
a 3 . The vacuum contribution to the decoherence func- 
tion is shown from Eq. (|20[) to be proportional to a. From 
the considerations above it follows that the emission of 
Bremsstrahlung photons doesn't seem to be relevant for 
the vacuum decoherence process. 

However let's observe that for short times (t <C fi _1 ), 
the decoherence factor T(t) of Eq. (|25|) and the mass vari- 
ation 8m(t) of Eq. ([44]) show both the same t and a de- 
pendence. This appears to suggest a connection between 
the decoherence process for small times (vacuum contri- 
bution) and the dressingprocess. In analogy to the case 
of the two level systems [2( , the link between dressing and 
vacuum induced decoherence could be attributed to the 
correlation that get established between each component 
p of the wave packet and the part of the dressing struc- 
ture of the transverse electromagnetic field associated to 
it. 

To verify this hypothesis, we shall analyze the evo- 
lution of the field associated to each component of the 
wave packet, during the initial phase of the decoherence 
process. 



A. Field structure dynamics 

In the analyses of decoherence the behaviour of the en- 
vironment is usually not investigated being the interest 
placed on the system evolution. In our case the environ- 
ment is the electromagnetic field and its behaviour during 
the decoherence process can be analyzed by performing 
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the trace of the total density matrix over the degrees of 
freedom of the particle. 

For calculation purposes we shall consider the initial 
wave packet of momentum width Ap as a sum of mo- 
mentum sharp wave packets of width Ap <C Ap. Each 
of these sharp wave packets is centered at a momentum 
p and it has in configuration space a width Af m H/ Ap 
taken less than c/VL so that the dipole approximation can 
be yet used. To describe the development of the field cor- 
related to one of these sharp wave packets centered at p 
we start from a totally decoupled initial condition. The 
field is taken in its vacuum state and the charged parti- 
cle is described by a sharp wave packet with momentum 
components peaked around p of the form ^/Wl J2 P $p,p\p) > 
where N e is a normalization factor and dp p indicates a 
quasi delta centered on p of width Ap. 

The corresponding initial density matrix is 



W) = fc(0)®l{05j}>({0y}l 

= J2 \P)NS^ P S^ P ( P '\ ® \{0- kJ })({0 kJ }\ 



(58) 



p,p> 



We shall consider the representation of p(t) in a coher- 
ent basis. Indicating with \\k,j) a coherent state of the 
mode {k, j} of amplitude A, the reduced density matrix 
elements of the field in this basis, with the initial condi- 
tion of Eq. (|58|) , are given by 



that compose the wave packet and the corresponding as- 
sociated transverse photons that leads to vacuum deco- 
herence in our system. 

To confirm the possibility of associating a number of 
photons to the various momentum components of a given 
wave packet, we could choose as initial state a sum of 
two sharp wave packets of width Ap peaked around two 
different momenta. In this case it is easy to show that 
the average number of photons surrounding the particle 
can be written as a sum of two terms relative to the two 
sharp wave packets composing the initial state. 

The energy associated to the field structure that builds 
up around the particle is responsible together with the 
interaction energy of the mass variation 5m computed in 
Eq. The average energy associated to the cloud of 

photons, obtained in Appendix ICl by (|C18|) . is equal to 



8a hvi n 2 t 



2+2 



p- 



3-7T moc 2 1 + Q 2 t 2 2mo 



(62) 



(Ep) can be written, using Eq. (|44|) for 5m, as (Ep) = 
— (p 2 /2m ) x (Spm/mo) with 5pm = —25m. Therefore 
(Ef) reflects on one side the build up of correlations with 
momenta and on the other side contributes to the mass 
variation. This explains the analogous time behaviour of 
5m(t) and T vac (t). 



= (*k,j\tvs{U(t)ps(0) ® \{0- kJ }) 

{{OkjKU- 1 ®}^,). (59) 

The explicit calculation, reported in Appendix [C] by 
Eq. (|C7[) , gives for the reduced density matrix of the field 



|A fcJ f \K',f\ 2 IC-I 5 



(60) 



with . defined in Eq. (|CJ2|) . Because of our choice of 
sharp wave packets in p space, the density matrix of 
Eq. (|60[) retains only a dependence on p. 

Eq. (|60"1) allows to get the average number of photons 
(hp) that can be associated to each sharp wave packet of 
width Ap and centered at the momentum p of the total 
wave packet. The calculation, performed in Appendix [Cl 
by (|C17[> . leads to 



2a 



P 



3ir m\(? 



In (i + fn 



2 ,2\ 



(61) 



The time dependence of the average number of photons 
of Eq. (fBTj) is, apart a factor 2, equal to that of the vac- 
uum contribution to the decohcrcncc function (|20p . This 
result appears to give a strong indication that it is just 
the buildup of correlations among the various momenta 



V. SUMMARY AND CONCLUSIONS 

We have considered a free charged particle interacting 
with a bath consisting of an electromagnetic field at tem- 
perature T. We have analyzed the decoherence on the 
charged particle wave packet induced by the interaction 
through the investigation of the off diagonal elements of 
the particle reduced density matrix. The interaction has 
been taken in the minimal coupling form and the particle 
is described by a wave packet of width Ar. The effect 
of all the modes of wavelength larger than Ar can be 
taken into account within the dipole approximation. The 
dipole approximation and the neglecting of the quadratic 
potential term reduces the coupling to a linear form and 
this in turn allows an exact treatment of the dynamics of 
the system. 

Our analysis has been conducted in the context of non 
relativistic QED which is in the spirit of modern quantum 
field theory an effective low energy theory with the cut 
off frequency £1 parameterizing the physics due to the 
higher frequencies [30l ]. For this reason our final results 
must show a dependence on Q, that is however as usual 
weak (logarithmic), as for example in the case of non 
relativistic expression for the Lamb shift. 

The analysis of the decoherence process has been con- 
ducted both in the momentum and configuration space 
and it has been possible to separate both the vacuum and 
the thermal contribution to decohcrcncc. 

In momentum space decoherence among different mo- 
mentum components occurs without population decay, 
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therefore decoherence occurs in its purest form that is 
without dissipation. This is reflected by the fact that 
the width Ap(t) of the wave packet remains constant in 
time while the coherence length l p (t) decreases in time, 
in particular as 1/Vt for large t. 



In configuration space again both vacuum and thermal 
contribution appear in the decay of the off diagonal el- 
ements of the reduced density matrix similarly to what 
occurs in the momentum space. However in the char- 
acterization of the development of decoherence by the 
behaviour of the space width of the wave packet, Ar(t), 
and the coherence length, l r (t), it is necessary to con- 
sider that in these quantities two contributions appear, 
which are not present in the momentum space. The first 
is due to the free evolution of the wave packet and the 
second to the dressing process. The appearance of these 
contributions only in the configuration space is due to 
the fact that the Hamiltonian commutes with each mo- 
mentum component that then results to be a constant 
of the motion. In particular the dressing process, with 
the emission and absorption of virtual photons and the 
creation of a structure of transverse field around the par- 
ticle, doesn't modify the distribution of momenta of the 
wave packet while it modifies the spatial probability dis- 
tribution. We have determined the contribution of these 
physical effects to Ar(t) and l r (t). 



We have tried to determine the physical effect respon- 
sible for the part of decoherence independent from the 
temperature. The Bremsstrahlung photons emitted dur- 
ing the dressing have been shown not to be relevant for 
vacuum decoherence. The results obtained about the 
particle mass variation indicate that the vacuum contri- 
bution to decoherence is temporally linked to the dress- 
ing process. We have shown by the analysis of the field 
structure dynamics that the onset of time dependent cor- 
relations, induced by the interaction, between the mo- 
mentum components of the particle wave packet and the 
associated field structure, may be held responsible of vac- 
uum induced decoherence. In fact the average number of 
entangled photons with a given momentum has the same 
time dependence of the vacuum part of the decoherence 
function and moreover has the same dependence on the 
physical parameters of the system. 



The results obtained for our system on the develop- 
ment of induced decoherence depend on the fact that in 
the initial state considered there are not particle- field cor- 
relations. Previously it has been shown that decoherence 
evolution is influenced by the presence of initial partial 
correlation [3l|, [H, [H, [3J] . It appears of interest to ana- 
lyze in which way the results obtained in this paper are 
modified in the more realistic case in which partial cor- 
relations between the system and the environment are 
present since the beginning. 



APPENDIX A 



Phase factor 



In order to compute the phase factor £ (p, t) of Eq. ([5]) it 
is necessary to explicit the commutator of Hi at different 
times ([7]). Using the following relation satisfied by the 
polarization vectors 



E m n _ sr jU ju 

hi kj u ra,n ™m n, n i 



(Al) 



where k indicates the versor of k and m and n are gener- 
ical components, we obtain 



j2\P-^\ 2 =P 2 -( P -kf 



(A2) 



By using the explicit form of the coupling coefficients 
of Eq. and Eq. (|A2[) to compute the sum over the 
polarizations, the commutator of Hi at different times 
([7]) assumes the form 

[6 I {a),A I {s>)] = -2i^fi*& p efa[{w k -k-v )(8-s>)] 

p,k,j 

Ane 2 hM-^p 2 - (p ■ k) 2 



mnV 



-tr p sin [(cjfc - k ■ vq)(s - s')} 



(A3) 



The time integrations present in Eq. ([5]) give 



t ,t 

ds / ds'sin [(uj k - k ■ v )(s - s')} 6{s - s') 



1 



sin(wfc — k ■ Vo)t 



(cjk-k - v ) 



(A4) 



and joining Eqs. (jAT|) . (|Al|) and © we obtain for t) 



2ire 2 \ - 



p 2 - (p ■ fc) s 



mlVh '-r' Lo k {u k - k ■ v ) 



sin^fc — k ■ Vo)t 



(ui k -k- v ) 



(A5) 



By taking the continuum limit on the field modes, 
V / °°d 3 fc/(27r) 3 , Eq. becomes 



Hm 2 (2n) 2 



d 3 fc 



sin(o; — k ■ v Q )t 
(u> - k ■ v ) 



p 2 x (l-k x )+p 2 y (l-k y )+p 2 z (l-k z ) 
x f —2— . A6 
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By introducing the cut off factor cxp(— w/O), Eq. (|A6|) 
assumes the form 



f° dw / 




/ — exp 


n) 


Jo u x 



(A7) 



do 



- X) - sin [u;£(l - X)] 



(l-X) 2 
+ pl(l-k 2 y )+pi(l- 



where do = sin 9d9dtp is the infinitesimal solid angle and 
we have posed 



fc ■ Vq = uj— [sin^o sin^fc cosf^fe — tpo) 
c 



COS Oq COS 6/fcJ 

= luX , (A8) 



where 8q and ipo are the angles of the vector Vo and 
and (fik are the angles of the vector fc. 
Indicating with f(X) the X dependent part within the 
integrand in Eq. (|A7|) . for small values of i>o/c, this can 
be expanded with respect to X obtaining up to the first 
order in X 



u>t(l - X) - sin [ut(l - X)] 



{i-xy 

ps - sincji + /' '(0)X . (A9) 
Using this expansion and the following integrals 
8 



J do(l-fe 2 



for i = x,y,z . 



(A10) 



y do(l - fe- )X = fori = a;,y,z, (All) 
in Eq. (|A7[) , we obtain up to first order in vq/c 



£(p,t) 



2eV 



2ap 2 
37rmgC 2 







/ — exp 




/o <*> v 



(fit — arctanilt) , 



j (w< — sin wt) 

(A12) 



where a = e 2 /he is a dimensionless coupling constant. 



Decoherence function 



To obtain the explicit expression of T P:P (t) (fTS)) it is 
necessary to calculate the trace on the field 



k,j *) = tr F { ex P 



p,p .,p,: 



p,p 



J*) P p } ■ 
(A13) 

The operator exp ( aj. i r y p ' P — afejlf'? *) is the generator 



of the coherent states of amplitude . It has been 

in 



shown [13, H3] that Eq, IA13l can be put in the form 

1 



Xilkj ,7ft j *) = exp 



(A14) 



Using Eqs. (|A13|) and (|A14|) . the relation 1 



trp 



= coth 



(A15) 



2k B T / 

obtained in the case of thermal distribution for pp and 

1 P ,p',2 u i\ -,2 [1 - cos(w fe - fe • vo)*] 47re 2 
l7 "' =Kp-*>'^ ^K-fe., ) 2 mp^ ' 

(A16) 

derived from the position following Eq. (|15p and using 
Eqs. (|3]) and ©, the decoherence function (fT"5|) can be 
put in the form 



27re 2 
tirrikV 



E 



x coth 



[1 — COs(Wfc — fc ■ Vo)t] 

ui k (u> k - k ■ v ) 2 



(A17) 



Taking the continuum limit on the field modes ^ fe — > 

Vj °°d 3 k/(2irf, using £. [(p-p') ■ £ fe ,,] 2 = (p-p') 2 ~ 

[(p — p') • fc] , inserting the cut off factor exp(— uj/fl) 
and introducing the variable X defined in Eq. (|A8[) . we 
obtain 



— exp 



--)coth 



2k B T 



do 



1 - cos [u)t(l - X)} 



(1 



k y 



X) 2 

+ (p-p')10- 



( P -p')l(i-k x ) 

(Ais) 



where, as before, do is the infinitesimal solid angle. 
Indicating with g(X) the X dependent part within the 
integrand in Eq. (|A18|) . for small values of vq/c, this can 
be expanded with respect to X obtaining up to the first 
order in X 

1 - cos [u>t(l -X)] 



9{X) = 



{i-xy 



1 - cos uit + g'(0)X . 



(A19) 

Using this expansion and Eqs. (|A10j) and (jAlljl to com- 
pute the angular integral in Eq. (|A18|) . we obtain up to 
first order in vq j c 



T p,p (f) = 



2a (p-p') 



'\2 



3- 



,2„2 



du 



■ exp 



x (1 - cos cot) coth (Hu)/2k B T) . (A20) 

Before carrying out the frequency integral in Eq. (1A20|I . 
we separate T™' (t) in two parts, TP'P' (t) = TP,£(t) + 
r?; p (t): a temperature independent part due to vacuum 
fluctuations and a dependent one due to the thermal bath 
properties, which goes to zero for T — > 0. From Eq. (|A20[> 
we obtain the temperature independent contribution as 

,'\2 



2a (p — p ' 



3tt 



dui exp 



2a (p — p 



>Y2 



3tt 



2 9 



In \fl + VlH 2 , 



(1 — cos Ult) 



(A21) 
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and the thermal contribution as 

x [coth (fiw/2k B T) - 1] 
2a {p-p'f 1 t* 



ds / dxexp 
3tt mgc 2 (3 J J 

x [coth (x/2) - 1] sin(sx//3) . 



kgTx 



(A22) 



For fc B r < fift and introducing r F = h/nk B T « 2.43 
KT 12 s/T[A], we find 

x [coth (x/2) — 1] sin(sx//3) 



2a (p — p ' 



In 



sinh (t/TF) 



t/rp 



(A23) 



37r mgC 2 

where in the integration on a; we have used the formula 

13 

dx [coth (2/ 2) — 1] sin(sa;//3) = 7r coth(7rs//3) . 

(A24) 

Summing the vacuum contribution given by Eq. (|A21|) 
and the thermal by Eq. (|A23[) , we obtain for the decoher- 
ence function 



T P;P ( f ) 



2a {p-p'f 

37T TO§C 2 



hi 



t/TF 



(A25) 



where we have introduced Ar = \/3d. q = r — t*o of 
components (x,y,z) and we have eliminated the explicit 
time dependence of T and Using 

[ + °° Pk ( b 2 \ 

/ exp(— ax 2 + bx)dx = \ — cxp ( — ) [Re a > 0] , 
J-00 V a \4aJ 

(B4) 

for the integral in p x we obtain 



h = F(p' x )- 



^pii exp 



„x2 ^ 2« a (r-^+^)-d 3 
~Po a «i(r-i*+4) 



(B5) 



x exp 



P qcI x 



4fi 2 (r-i$ + f) fi3(r-«$ + f)_ 

where F{p' x ) contains the integral in j4 and is equal to 



F(jp' x ) 



dp' x exp 



jr-z<i> + §| 2 -r 2 



, / 2d 2 pl ix' iHTx + 2p*d 2 T 



h- 

xj2i 



h 2 H h 2 (T - i$ + &) 



^(r-^ + f) 



2r 



d 2 , d 4 
7^ + 7F 



<J>2 



(B6) 



cxp' 



2p%d 2 {2Y + |s - ift(T - i$ + |*)a;' + iftrx 



4fi 4 (|r-z$ + f | 2 -r 2 )(r- l( i> + f ) 



APPENDIX B 

Here we report the explicit computation of the spatial 
reduced density matrix that involves the double Fourier 
transform of the reduced density matrix in the momen- 
tum space: 



1 



(27Tfi) : 



d^dVPs P (t)exp 



-{p-r-p r ) 



(Bl) 



Using Eq. ^ for p p s ' p (t), with the initial wave packet 
form of Eq. (|30|) . we can easily decompose Eq. (|B1|) in 
equal components: 



P7 (*) 



N x 4 N y / y N z /, 



(27rfi) (27rfi) (2?rfi) 
where x,y, z are mute indices, N x = l/v^TrApx and 
2d 2 p% 2 



(B2) 



Jo; = exp - 



fi 2 



exp 



Jdp x dp' x x (B3) 
2d 2 pg ix 



- P 2 x [r- i -f + ^)+p x l2Tp' x , fi2 



x exp 



-p'i r + 2 $ 



fi 2 



, / 2d Pq ix' 



Substituting this result in Eq. (|B5[) . simplifying and ra- 
tionalizing where it occurs, and posing the adimensional 
quantity Z = 1 + 2Th 2 /d 2 + $ 2 H 4 /d 4 , we obtain after a 
lengthy calculation 



irH 2 

It. = 1= : exp 



Zrf 2 



pg 2 d 2 (2r - i2$ + £)(T + z$ + 



x exp 



x exp 



x cxp 



pg 2 (2r-z$+^) 2 (r + z$+§) | 2rx:c' 

(r 2 + z^)z 

(r + i$ + ^)a; 2 + (r-i$ 



€)x' 2 



AZd 4 /H 2 

: (2r + z$ + ^) P gx-(2r- t $ + ^) P gx' 

' Zd?/h 



AZd 4 /h 2 

(B7) 



After some passage, from I x we obtain p r ^ r (t) 
NI x I y I z /(2irH) 3 put in the form 



Ps (*) = „ , Kr7 3 exp 



(l + 2rfi 2 /d 2 )p -( (7 -9 / ) 



x exp 



x exp 



3[(g + 2 Po $ft) 2 + (g' + 2 Po <M) 2 
4(3d 2 Z) 

3rfi 2 (g ~q') 2 /d 2 + i3$fi 2 (g 2 - q' 2 )/d 2 
4(3d 2 Z) 



(B8) 
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The last can be put in a useful form to compute directly 
some quantities, as 



PT (*) 
x exp 

x exp 

x exp 



= s¥zl exp [ mpz 

.(l + 2Th 2 /d 2 )p Q -(q-q') 



(B9) 



hz 

{l + 2rti 2 /d 2 ){q-q') 2 + l{q + q^ 



A(3d 2 Z) 

-12$%> • (q + q')-i3<$>h 2 (q 2 -q' 2 )/d 2 ^ 



4(3d 2 Z) 



APPENDIX C 

We compute the trace on the subsystem in the mo- 
mentum basis, \p(t)). In the interaction picture we 
have \p(t)) = exp[itp 2 /2m n H]\p(0)), but because here 
we are not interested in the free evolution we will use 
\p(t)) = \p(0)) for the trace, thus neglecting the phase 
factor which isn't relevant for the following discussion. 

We rewrite the time evolution operator of Eq. © in 
the form 

u(t) = Yl cx p K(p, t)v P ] cx p °v - 

P.k,j 

(CI) 

where (3 k ■ = g k -au is given, using Eqs. ((4]) and ([9]), by 



Pk,- 



m V Vu)k fi{oJk — k-vo) 
Using Eq. (|UT|) in Eq. $5§§ we obtain 



(C2) 



Pf ' 3 ' Afc ' ,y (t) = (Afej|tr s 



n n 

Pi,ki,ji P2 : k 2 J2 



{exp \a Pl (&l uh pX u h ~ ^xP%\*jij\ CX P [*f(pii*)*Pi] 
x E \p)\{^j})({OhMP'\ N ^P S p',P 



p,p 

x exp \a p2 {-al 2d JH h + h 2 ,jj] 
x exp[-i£(p 2 ,t)o-p 3 ]}} lAfc/^,) . 

Using the cyclicity of the trace and 
exp (T p 



,p 2 * 



(C3) 



= \p) cxp (alt ./?^ - a fcJ /3£*.) I« fcj > = |P>|^> , (C4) 
where the amplitude of the coherent state \0f. j) depends 



on the momentum component p, Eq. (|C3[) becomes 



Pf 



(*) = / dVX^-PW-^ 



>< n n (^i^x^jj^j/) • (C5) 

fcl Jl &2 J2 

Taking into account the explicit form of the scalar prod- 
uct between coherent states 



k',j' 



I-^jI 2 l^fcjl 2 , x * r p 
— — + A fc,iP feJ 



(C6) 



Eq. (|C5[) becomes 



Pf 



(t) = N £ exp 



\Xkj\ 2 \K;y? l/?Lf 



(C7) 



where the integral over the momenta has lead to the pres- 
ence of p in the (i v k ■ . 

Now, we calculate the average of the operator aj. -&k,j 
using the trace in the coherent states basis, that for a 
generic operator A has the form (2^1 



TrA = — I d(RcA)d(ImA)(AU|A) 

7T 



(C8) 



Using a^ j&k,j as operator in Eq. (|C8[) we obtain (omitting 
the pedici k and j in A) 

= Tr F (alj&kjPF) = Tr F (&\ >:j PF&k,j) 

- Ti F ( PF ) = ^ / d(RcA) d(ImA) (|A| 2 - l) p x / , 

(C9) 

where we have used the commutation rules satisfied by 
the operators h\ j and a,k,j, and the action of these oper- 
ators on the coherent sates 



a fej |A)=A|A) and (A| a£ , = A*(A| 



(CIO) 



Substituting the diagonal elements of Eq. (|C7[) in 
Eq. (|C9|) . we obtain 

=\ j d (RcA) d(ImA) (|A| 2 - 1) (Cll) 



x exp 



-|A| 2 -|/3LI 2 + ^L + ^ 



Posing RcA = Ai and ImA = A2, Eq. (|C11[) can be put in 
the form 



exp(-|/3f,,| 2 ) 



dAi dA 2 (A 2 + Al - 1) x 



(C12) 



exp 



-A 2 + 2Ai {Ref3 p k j 



exp 



-A? + 2A 2 
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The integrals involved in Eq. (|C12[) are of the Gaussian 
type (|E4|) and 



dx x exp( 



-ax 



bx) 



4a 5/2 



■ exp 



4a 2 

* (C13) 

Using Eqs. (|C13[) and (|B4)) in Eq. (jCllj) we obtain easily 



Using the expansion of Eq. (|A19[ ). and Eqs. (|A10[) and 
(|Alip in Eq. (|C16[) . we obtain that the angular integral 
is equal to (1 — cosujt)8TTp 2 /3, while the resulting inte- 
gral over frequencies gives In (l + VL 2 t 2 ) /2. The average 
number of photons at time t that compose the cloud as- 
sociated to momentum p is then equal to 



Aire 2 (p ■ Ekj) 
= Vm 2 hu;l(l-X) 2 
x {1 - cos [w k t(l - X)]} , (C14) 



with X defined in Eq. l|A8p . Eq. (|C14|) represents the 
average number of photons, of the mode of the field rep- 
resented by {k, j}, that compose the cloud associated to 
the momentum p. To obtain the trend of the total num- 
ber of photons, then, we must sum over the polarizations 
and over k. Using Eq. (|A2[) to perform the sum over j 
and Eq. (|A8[) we obtain 



£<4 



fc,7" afe J/ 



(p-k) 



ine 2 



l-cosM(l-X)] 

x {T~xJ 2 ■ (C15) 

Performing directly the limit to continuum on the field 
modes and inserting the usual cut off factor exp(— lo/VL), 
the sum over the fc assumes the form 



V 



(2tt) 3 

,,2 



/d 3 fcex P 



2tt 2 



m 2 hc 3 ,/Q 



duj 



■ exp 



d<7 



1 - cos [ut(l - X)] 



(1-X) 2 



f x (l-kl)+pl{l-kl 



Pin 



(C16) 



2q 



3ir uiqC 2 



\n(i + n 2 t 2 ) . 



(C17) 



The average energy associated to this cloud of photons 
is obtained multiplying the integrand of the final integral 
over the frequencies (|C16jl for Huj and using considera- 
tions similar to those leading from Eq. (|C16j) to (|C17I) : 



(Ef) 



V 



(2tt) 3 
4a Hp 2 



d k exp 



a/c. 



Sir ttIq 
8a hVL 



%-rr \ dwexp (-^) (1 - COSLUt) 
o c Jo v " y 



2+2 



P 



3-7T moc 2 1 + Q 2 t 2 2mo 



(C18) 
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